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ABSTRACT 

In this paper, we present a new numerical formulation of solving the Boundary 
Integral Equations reformulated from the Helmholtz equation. The boundaries of the 
problems are assumed to be smooth closed contours. The solution on the boundary is 
treated as a periodic function, which is in turn approximated by a truncated Fourier 
series. A Fourier collocation method is followed in which the boundary integral equation 
is transformed into a system of algebraic equations. It is shown that in order to achieve 
spectral accuracy for the numerical formulation, the non-smoothness of the integral 
kernels, associated with the Helmholtz equation, must be carefully removed. The 
emphasis of the paper is on investigating the essential elements of removing the non- 
smoothness of the integral kernels in the spectral implementation. The present method 
is robust for a general boundary contour. Aspects of efficient implementation of the 
method using FFT are also discussed. A numerical example of wave scattering is given 
in which the exponential accuracy of the present numerical method is demonstrated. 


‘This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NASI- 19480 while the author was in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23681. 
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1. INTRODUCTION 


The Boundary Integral Equation Method is a powerful tool for solving certain 
boundary value problems. It is particularly attractive in developing numerical methods 
since, when applied, it reduces the dimension of the problem and often transforms a 
problem in an infinite domain to integrals on the finite boundary in which the far field 
radiation condition is satisfied automatically. 

Numerical methods for the Boundary Integral Equations have been developed 
predominantly based on the Boundary Element Method (BEM). In this method the 
physical boundary is divided into finite elements and integrations over each element are 
approximated by numerical quadratures. In this way, the integral equation is converted 
into a system of algebraic equations. BEM has gone through a rapid advancement in 
recent years. Its applications to the Helmholtz equation are discussed in reference [2] 
and the references cited therein. Other formulations of solving the Boundary Integral 
Equations in wave scattering and propagation have also been proposed in the past, 
including, for example, the T-matrix Method[3,4,5j. 

In this paper, we present a spectral method formulation for solving the boundary 
integral equations reformulated from the Helmholtz equation. For the Fourier approx- 
imations used in this paper, we assume that the boundary of the problem is a smooth 
closed contour. The functions appearing in the 2-D boundary integral equation will 
be treated as periodic functions, which are in turn approximated with high accuracy 
using truncated Fourier series. The boundary integral equation is then transformed 
into a system of linear algebraic equations. The spectral methods have been known to 
have extremely feist convergence rates, faster than any finite power of 1 /N when the 
solution is infinitely smooth[6,7] (where N is the number of collocation points). The 
present numerical formulation will be seen to have spectral accuracy. 

It is known that, although any periodic function can be approximated by a trun- 
cated Fourier series, the rate of convergence depends on its smoothness. Unfortunately, 
the integral kernels for the Helmholtz equation are not smooth. In particular, we note 
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that the 2-D Green’s function, appearing in the integral equations, possesses a loga- 
rithmic singularity. Furthermore, the normal derivative of the Green’s function also 
contains a term involving the logarithmic function. It will be seen that it is critical to 
remove the non-smoothness of the integral kernel in order to achieve fast convergence 
in the spectral formulation. In this paper the non-smoothness of the integral kernels 
is subtracted out by using a logarithmic periodic function whose Fourier expansion is 
known. 

Fourier approximation methods for the boundary integral equations in wave scat- 
tering have been proposed in the past for simple geometries. A “fast numerical method” 
has been formulated by Bojarski[8] for wave scattering by a circular cylinder. It was 
pointed out that the boundary integral equation for scattering by a hard circular cylin- 
der can be solved easily and efficiently in the Fourier spectrum domain of the solution. 
Due to the special geometry considered, a simple relation for the Fourier coefficients 
of the solution and those of the boundary condition was derived. Recently a similar 
numerical approach has been applied by Schuster[9] for a wave transmission problem 
of concentric cylinders. However, in these works the singularities and non-smoothness 
of the integral kernels were not removed. Consequently the convergence rates of these 
methods were not exponential. 

In section 2, a formulation of the Boundary Integral Equation for the Helmholtz 
equation is given, followed by a discretization using a Fourier collocation method. 
Essential elements of achieving spectral accuracy are investigated. In section 3, a 
numerical example is given in which the spectral rate of convergence is demonstrated. 
Section 4 contains the conclusions. 


2. FORMULATIONS 

2.1. Boundary Integral Equations 

Consider wave propagation in an interior or exterior domain with a smooth closed 
boundary I\ Figure 1. The wave equation for a scalar function <f> with assumed time 
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dependency of e ,wt is reduced to the Helmholtz equation 


vV + «V = o 


(i) 


d 2 


where k is the wave number and V 2 is the 2-D Laplace operator, V 2 — + 


The boundary condition considered in this paper is one of the following types : 


dy 2 


Dirichlet : 4 >(r) = a(r') On T 

Neumann : ^(r) = b(r) On T 

on 

Robin( Impedance) : a + /? -^(r) = c(r) On T 

The Helmholtz equation ( 1 ) can be reformulated into a Boundary Integral Equa- 
tion in various ways[ 10 ,ll]. Our purpose in this paper is to demonstrate the numerical 
method using spectral approximations and the essential elements of achieving spec- 
tral accuracy of the numerical solution. A direct formulation of the Boundary Integral 
Equation employing the Green’s Function will be used here which leads to an integral 
equation involving <j) and its normal derivative on the boundary [2]: 



where 7— denotes the normal derivative, with the direction of n being outward from the 
on 

domain of wave propagation, and r*p denotes any point on the boundary. The Green’s 
function G(r, r'), whose form will be given later, satisfies the following equation 


V 2 G + K 2 G = -6(r-r') 


( 3 ) 


Once the values of 4 > and 77— on the boundary are found, the solution at any point, 

on 

r, inside the domain of wave propagation can be obtained as [2] 



3 



Without loss of generality, let us now assume the boundary curve T be parame- 


terized as 


f r = r(0), 0 < 0 < 2tt 


where 0 is a non-dimensional parameter for the boundary contour, not necessary the 
arclength or some angle. In addition, the curve is supposed to be simple, that is 

For simplicity, we further assume that r(9) is infinitely differentiable. Then the bound- 
ary integral equation (2) can be written as 

!*«-)+ jfwffr*-) ||| *- jf H| -» w 

dr 

in which the line element dT = — d9. For clarity, the dependency on 6 has been 

d0 

expressed explicitly in (5). 

For the 2-D Helmholtz equation, the out-going Green’s function and its normal 
derivative are[12] : 

G(6,0 l )= t -H i o '\ K R) (6) 


G(6,0') = 1 H ( q 1, (kR ) 

(Q o>\ _ * K ‘ n 


in which we have used the notations that R = r(6) — r*(d / ) and R = |JJ|. Here and 
H™ are the first kind Hankel functions of order zero and one, respectively. 


2.2. Spectral Approximations 

Q / 

For a closed boundary T, <f>(6) and — (0) are periodic functions of 0. for 0 < 0 < 

d* . 

27r. Let <f> and — - be approximated by truncated Fourier series as : 
on 

Nf 2-1 

*"(«) = Y ( 8 ) 

n= — N/2 
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( 9 ) 


*/2" 1 

Vy ( a\ V' J, „«»»* 

»r (9) = 2^ *- e 

n=-N / 2 

The particular form of the truncated Fourier series has been taken in ( 8 ) and 
(9) for the convenience of using discrete Fast Fourier Transform (FFT) programs[7]. 
Substituting ( 8 ) and (9) into (5), the boundary integral equation becomes the following 


N/ 2-1 


AT/2-1 


iV/2-1 


*t/r» \r /o *2 U AT /O ^ u 


n— — N/2 


n=-N/2 


n=-N/2 


We note that the two integrals in (10) are actually the Fourier coefficients of 
dr OG dr 

G(0, 9’) — and (9, 8') — respectively. Our aim is to evaluate the integrals by 
d9 on ad 

the Fast Fourier Transforms with spectral accuracy. 

dr* dG dr 

It is clear that both G(9,0 f ) — and (0, 6') — are periodic functions of 6 

du on du 

and 9 f . They are also infinitely differentiable except at 9 = & where R = 0 in ( 6 ) 

and (7). Although any periodic function can be approximated in a truncated Fourier 

series, the rate of convergence depends on the smoothness of the function. From this 

consideration, we note that, first, G(9 y 0 i ) has a logarithmic singularity at R = 0 , due 

dG 

to the Hankel function of order zero in ( 6 ). Second, although - 5 — (0, 8 f ) can be shown 

on 

to be a finite function, it is not infinitely smooth. It is easy to show that the function 
R has a discontinuous derivative at 0 = 9*. In particular, for \9 — 0 ; | small, we have 

*-W)-iWI = |f •(»-»-) + ! ^ ■(»-•’)* + 

dr 1 dP r 

(11) 

Consequently any term with an odd power or logarithmic function of R will not be 
infinitely smooth and has to be treated in order to achieve spectral accuracy in the 
Fourier approximation. 
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To study the singularity in G(9, 9'), we note that 


G(9,9') = = 1 -(J 0 (kR) + iY 0 ( K R)) 

in which Jo and Yo axe the zeroth order Bessel functions of the first and second kind, 
respectively. Using the asymptotic series for a small argument[13], 




Jo(kR) 


Yo(kR) = - In J 0 (kR) + ^Jo(kR) + ^ 

7T \ Z J 7T 2 


we have, for \9 — 9'\ small, 


G(9,9') = l -(J 0 (KR)+iY 0 (KR)) 


= — — In Jo(kR) + ( smooth terms 


in which “smooth terms” represents a convergent power series containing only the 
even powers of JZ, and Jo is the regular Bessel function of the zeroth order [13], also a 
power series of R 2 . To remove the logarithmic singularity but preserve the periodicity, 

0 _ 0t 

consider the periodic function In 2sin( — - — ) , which has a logarithmic singularity at 
9 = 9'. Its Fourier series can be found exactly [12] : 

i . , 9 — 9' cos n(9 — 9') 

In 2 sin( — ) = - £ L i (12) 


To subtract out the singularity in G(9,9'), let 

<5(0, S') = i-ln 2sin(— ^— ) Jo(-cfl) 

Then the Green’s function is found as 


G(9,9') = 0(9,9') - i-ln 2sin(i^) J 0 (kR ) (13) 

Using (11), it is easy to show that 0(9,9') is periodic but finite for all values of 
the arguments. Furthermore both 0(9,9') and Jq(kR) are now infinitely differentiable 
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functions. Thus the Fourier coefficients of G(0, O') — can be computed according 

dO 

to (13) with spectral accuracy. For the term involving the logarithmic function (12), 
convolution sums will be used. 

We now study the singularities in the normal derivative of the the Green’s function 

dG 

- 5 —. The asymptotic expression of h[ X \kR) for kR small has the following form[13] 

on 

H^inR) = Ji(kR) + iYx(KR) 

= Ji(kR) + i ( + — In ( Ji(kR) + odd powers of kR ) 

\ 7 tkR 7T \ 2 / / 

2 i 2 i . ( k R\ . 

= — H In — Ji{kR) + odd powers of kR 

7r kR 7 r V 2 / 

in which we have used the fact that the Bessel function of order unity, Ji, is a power 
series of odd power of kR, 

r / r.\ K R (kR) 3 

«' R )=T~ 16 +- 

In addition, it can be shown that R • n = O((0 — O') 2 ). Thus it follows that 

dG .. ... ik ( 2 i 2 i , ( kR\ , , . R • n 

= + 7 ln( T )./,(Kfl)+ odd Powersoft) — 

k , / /ci2\ _ . n .R- n , . 

= — In I — ) Ji(kR ) — — — |- ( smooth terms ) 

27T \ 2 / XL 


It is seen that, although is a finite function, it does not have smooth derivatives due 

On 

to the logarithmic function appeared in the first term shown above. For this reason, 
its Fourier approximation will converge only at a finite rate of 1/iV 3 . To smooth out 
the function for computations by FFT, let 




0 - 0 ' 


H{0,0') = -jHY\kRY-^-^- In 2 sin( — ) *(«*) 


Then 


?|W) = «(«.«■) 1” 2sin(^— -) 
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- r a 

The functions H(8,8') and J\(kR) — — — in (14) are now periodic and infinitely 

R 


dG 

differentiable. The Fourier coefficients of -r— (6, 6') 

on 

(14). 


dr 


de 


will be computed according to 


Following the above analysis, we then express the boundary integral equation (10) 
as follows : 


1 N/2-1 . 2n 

o E #•«"*'+ E M / 


N/ 2-1 


n=-N/2 


n=-N/2 


dr 


do 


d$ 


+ 


— [ 2 * , 
2tt J 0 


e in6 In 


n . ,0-0\ r , ^R-n 
2sm (^— ) 


dr 


d6 


d6 


N / 2 ~' / />2ir 

= E M / c ‘ ne 

n=-N/2 ' J ° 


G(M') 


dr 


d8 


d8 


_ _L [ 2 * 

2 7T J 0 


e‘" tf In 


„ • ,*-0\ 
2sin (— 2~ ) 


Jo(kR) 


dr 


d8 


<») 


(15) 

Now all the integrals in (15) axe in a form which can be evaluated numerically 
with spectral accuracy. 


2.8. Discretization 

In what follows, a spectral collocation approach will be taken in deriving the 
algebraic system for the boundary integral equation (15). Let us introduce two sets of 
discretization points (Figure 1) 

8j = 2ttj/JV, j = 0, 1, ...,N — 1 

ff, = 2*0' + 1)/1V, ;' = 0,l,...,Ar-l 

The reason for using two sets of discretization points, as will be clear later, is that it 
avoids the direct evaluations of G(8, 8') and H(8,8') at points where 8 = 8'. Although 
both functions are finite there, their limits are geometry dependent. The current dis- 
cretization is robust. 
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We then require that the boundary integral equation (15) be satisfied at collocation 
points O' = d'y, j' = 0,1, ...,N — 1 . 

For convenience of discussion, we denote the following Fourier series approxima- 
tions at 0'y for j' = 0, 1, N — 1, 


Ji(kR) 


G(9,9'y) 

3(0, 9'y) 

[Jo(kR)]v" 

R • ft 


dr 


d9 

dr 


dO 

dr 


R 




dO 

dr 


d9 


N/ 2-1 

E —I TtO 

gj'nC 

n=-N/2 

N/2-1 

= £ h i'nt~ in9 

n=-N/2 

N/2-l 

= Ys Pj'n ^~ ,n9 

n--N/2 
N/ 2-1 

_ V „ -in# 

— 2^ gync 

n=-N/ 2 


(16a) 

(166) 

(16c) 

(16d) 


The coefficients of the expansions are calculated by FFT (backward in the usual 
sense) as follows: 

N-i 


9 ? « = 


>=o 

N - 1 


dr 

dd 


(«>) 


k = j E 


j= 0 

N - 1 


dv 


Pi ' n = jy £ [*(«*)],;, 


i=o 

iV-1 


dr 




i=o 


Ji(kR) 


R ■ ft 
R 




dr 

Te 


(16a') 

(166') 

(16c') 

(16d') 


In addition, we denote (12) as 
In 


2sin (— 2~ ) 


— ^ a„ 




(16e) 


1 


where a 0 = 0 and a n = j — - for n #0. 


n 
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It is immediately seen that the two integrals involving G and H in (15) equal to 
2^9 j*n &nd 2 tt hj* n respectively. The other two integrals are obtainable by convolution 
sums. Specifically, using the definitions given in (16), we have 

AT/2-1 

Y (17a) 

m=-N/ 2 

N/ 2-1 

<*>= y i/-"-..* 1 - 

m= — N/2 

By substituting (16) and (17) into (15), the resulting algebraic system is : 

1 Ar/2 ~ 1 , , N/2 ~ l A N/2-1 

2 S ^e n6 i' + <^n (2nhj> n + KUj' n ) = ^2 'Pn(2irg j '„-v j . n ) (18) 

n=-N/2 n=-N/2 n=-N/2 

for j' = 0, 1,2, N — 1. 

Equation (18) is easily cast into a matrix form 

A$ = B4> ( 19 ) 

where $ and are the vectors containing 0 n and r/>„, respectively. For Dirichlet 
boundary condition, 4* is solved from (19) with $ obtained from the boundary condition 
and vice versa for the Neumann boundary conditions. The Robin type impedance 
boundary condition can be treated similarly. 

2-4 Evaluation of convolution sums 

The convolution sums in (17) require 0(N ) multiplications for each uy n and vy n . 
Thus the total operations of obtaining the convolution sums are of order 0( N 3 ). This 
cost can be reduced considerably to 0(7V 2 log 2 N ) by the use of a pseudospectral trans- 
formation method with de-aliasing techniques[6,7]. For completeness, the evaluation 
of (17) with a “padding” de- aliasing technique is given below. 

Let M > 3 N and 

(j = 2ir j/M , j = 0,1,2, ...M-l 



Uj ' n = hJ 0 e ‘ n ^|2sin(-2^)| Jo(kR) ||| dS = 

—if •“ 4 *< i r fc » «■*>¥ S 
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Compute the following using FFT for j = 0, 1, 2, ...M — 1 : 


Af/2-1 

Aj'j = ^2 e~ xmii 

m~ — M/2 

M/2-l 

P,;= Y. 

m= — Af/2 


where 


«m = 


Pj'm = < 


a m —A 7 ’ < m < iV — 1 

0 other 

f Pj'm —N/2 <rn < N/2 - 1 
0 other 


and form the product 

Uj/j = AjijPjij 

Then the convolution sum Uj> n is the (backward) FFT of Uyj : 

1 M-l 

“>■» = T? E 


j=0 


for —N/2 < n < iV/2 — 1. The convolution sum can be obtained identically. 


3. NUMERICAL EXAMPLE : SCATTERING OF A PLANE WAVE 
BY AN ELLIPTIC CYLINDER 

The numerical method described in the previous section will now be applied to 
the problem of a plane wave scattering by an elliptic cylinder. Exact solutions of wave 
scattering by an elliptic cylinder can be easily obtained in infinite series[12]. However, 
numerical evaluation of the exact infinite series is not easily obtainable due to the 
complexity of the Mathieu functions involved. Our purpose here is to demonstrate 
the exponential rate of convergence of the spectral method formulated in the previous 
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section. The numerical results will also be compared with some asymptotic values 
available in the literature[14]. 

S.l Spectral rate of convergence 

Let the plane incident wave having an incident angle a about the x axis be 

^i/c(x cos Qf-f y sin a) 

The scattered wave, <f>, satisfies the Helmholtz equation (1). The boundary conditions 

considered here are the Dirichlet (soft) type <j> = —<f>i or the Neumann (hard) type 
d<f> d<j>i 

dn dn 

A simple parameterization of the elliptic cylinder is given by 

T : (x,y) = (a cos# , 6 sin#) , 0 < # < 2n (20) 

where a and b denote the major and minor axis of the ellipse. The normal vector on T 

is n = —(b cos #, a sin #)/ y/tf cos 2 # + a 2 sin 2 #. 

Numerical results for the Dirichlet boundary condition will be presented first. For 

this calculation the parameters have been taken to be a = 2, b = 1 and k = 2n. In 

table I, the values of the solution at three selected points on the boundary are given 

as the number of collocation points increases. Listed are the values of — — at # = 0°, 

on 

90° and 180° on the boundary for the incident angle o = 0. Since no exact value is 
available, the numerical solutions are compared with the results obtained for N = 256. 
Clearly exponential rate of convergence is observed. 

We point out that although the parameterization of the ellipse given in (20) is 
a convenient one, it may not the be the best. For example, the following form of 
parameterization, which has more uniformly distributed collocation points than (20), 

T: (x,y) = cos (# + 0.05 sin(2#)) , 6 sin (# + 0.05sin(2#))^ , 0 < # < 2-n (21) 

yields the values shown in Table II. A somewhat better accuracy for small N is observed. 


12 



Q1 

In Figure 2 the solution, — , as a function of 8 on the surface of the elliptic cylinder 

on 

has been plotted. Shown are the results for the number of collocation points being 32 
and 256. 

8.2 Far field directivity 

Fax field scattered intensities, computed as r^ 2 , have also been calculated. The 
values on an exterior point are obtained through the boundary integral as 



The Green’s function for points lying outside of the boundary does not have any sin- 
gularity. Thus (22) is evaluated using FFT directly. The directivities of the scattered 
intensity for wave incident angle a = 0°, 30°, 60° and 90° are shown in Figure 3. For 
this calculation, a = 2, b — 1 and k = 2n. 

Barakat[14] gives some asymptotic values of long wave scattering of an elliptic 
cylinder at far field for normal incidence. A comparison of the current numerical 
results with the asymptotic values is presented in Tables III and IV. For these cases, 
a = cosh /io , b = sinh fiQ and k = 0.2. Here /zo corresponds to a value in the elliptic 
coordinates. The results of the far field scattered intensity for fio — 1.0, 0.6 and 0.2 
are presented in Table III and IV for the Dirichlet and Neumann boundary conditions, 
respectively. For the present long wave scattering, N = 64 has been used for all the 
calculations. The numerical computation agrees with the asymptotic estimation. 

Finally we note that although the exterior scattering problem is uniquely deter- 
mined, the direct formulation of the Boundary Integral Equation used here would not 
yield a unique solution when the wave number k coincides with the eigenvalues of a 
correspondent interior homogeneous problem. On the other hand, this numerical diffi- 
culty is well understood and remedies are readily available[10,ll]. For example, unique 
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solution can be obtained in a combination of single- and double-layer formulation of the 
Boundary Integral Equation[15]. The modifications of the present spectral implemen- 
tation to other formulations of the Boundary Integral Equation are straightforward. 

4. CONCLUSIONS 

In this paper, a spectral method of solving the Boundary Integral Equation is 
presented. It is shown that the integral kernels for the Helmholtz equation contain 
singular terms that have to be removed to achieve the spectral accuracy. Detailed 
numerical implementation of a Fourier collocation formulation has been given. The 
non-smoothness of the integration kernels is subtracted out by using a logarithmic 
function whose Fourier expansion is known. The numerical formulation presented here 
preserves the spectral accuracy and yields an exponential rate of convergence. 

Compared to the Boundary Element approaches, the Spectral Boundary Integral 
Equation Method presented in this paper would yield matrices of much smaller size 
since the latter requires far fewer points to achieve the same accuracy. This, of course, 
reduces the complexity and the cost of solving the resultant algebraic system. Since 
both methods will result in dense matrices, it appears that the spectral formulation is 
more advantageous for problems with smooth geometries. 
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Values of 


d<j> 

dn 


TABLE I 

for the Dirichlet boundary condition at selected points on the boundary. 

The boundary is given by (20). 
a = 2, b = 1, k = 27t, incident angle a = 0°. 


N 

II 

O 

o 

O 

O 

II 

6 = 180° 

Relative Error 

32 

6.178009567 

2.376671291 

7.134225119 

0.04 

48 

6.379285016 

2.175324633 

6.914464114 

0.2xl0 -4 

64 

6.379334010 

2.175259123 

6.914387795 

O.lxlO -7 

80 

6.379333962 

2.175259114 

6.914387773 

10~ 9 

96 

6.379333961 

2.175259114 

6.914387772 

lO" 11 

256 

6.379333961 

2.175259114 

6.914387772 



Values of 


d<f> 

dn 


TABLE II 

for the Dirichlet boundary condition at selected points on the boundary. 

The boundary is given by (21). 
a = 2, b = 1, k = 27r, incident angle a - 0°. 


N 

6 = 0° 

9 = 90° 

6 = 180° 

Relative Error 

32 

6.402635098 

2.165098187 

6.899520927 

0.005 

48 

6.379355527 

2.175263759 

6.914393819 

0.3xl0 -5 

64 

6.379334276 

2.175259159 

6.914387893 

0.5xl0 -7 

80 

6.379333966 

2.175259114 

6.914387774 

0.5xl0 -9 

96 

6.379333961 

2.175259114 

6.914387772 

10~ n 

256 

6.379333961 

2.175259114 

6.914387772 
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TABLE III 

Values of rtj> 2 at fax field for the selected angles. Dirichlet type 
soft wall boundary conditions are applied. 
a = cosh fio, b = sinh^o? « = 0.2 and incident angle a = 90°. 


Fax Field Angle 

Numerical 

Asymptotic[14] 

Ho = 1.0, 9 = 0°,180° 

1.769057417 

1.774 

fi 0 = 1.0, 9 = 90° 

2.137413140 

2.113 

O 

II 

O 

Qb 

II 

to 

-a 

o 

o 

1.538098511 

1.491 

Ho = 0.6, 9 = 0°, 180° 

1.364415306 

1.367 

Ho = 0.6, 6 — 90° 

1.506576730 

1.484 

Ho = 0.6, 9 = 270° 

1.287260683 

1.258 

Ho = 0.2, 9 = 0°, 180° 

1.062775602 

1.064 

H o = 0.2, 9 = 90° 

1.107062467 

1.087 

Ho = 0.2, 9 = 270° 

1.062022484 

1.041 
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TABLE IV 

Values of r<f> 2 at far field for the selected angles. Neumann type 
hard wall boundary conditions are applied. 
a = cosh/io, b = sinh^xo, « = 0.2 and incident angle a = 90°. 


Far Field Angle 

Numerical 

Asymptotic[14] 

Ho = 1.0, 9 = 0°, 180° 

0.009175011010 

0.00884 

Ho = 1.0, 9 = 90° 

0.02258893859 

0.0216 

Ho = 1.0, 9 = 270° 

0.1110188572 

0.1113 

o 

O 
00 
T— I 

o 

O 

II 

<q 

o 

II 

o 

3. 

0.001676516587 

0.00165 

Ho = 0.6, 9 = 90° 

0.007209616977 

0.0071 

Ho = 0.6, 9 = 270° 

0.02745548914 

0.0273 

Ho = 0.2, 9 = 0°, 180° 

0.0001288604229 

0.00013 

Ho = 0.2, 9 = 90° 

0.003710427872 

0.0037 

Ho = 0.2, 9 = 270° 

0.006989684179 

0.0069 
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<f 


Figure 1. Schematic of discretization of the boundary. 6 is the parameterizing variable 
for the boundary contour. n e and n, indicate the direction of the normal vector for an 
exterior and interior problem respectively. 
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d<f> 

Figure 2. Solution — on the surface of the ellipse for the Dirichlet boundary condition. 
on 

a = 2, b = 1 and k = 2n. Plane wave incident angle a = 0°. The parameterization of 


the ellipse is that of (21). 
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Figure 3. Far field directivities of the scattered wave for indicated incident angles. 
a = 2, b = 1 and k = 27r, with Dirichlet boundary condition. 
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